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We investigate the cosmological dynamics of non-minimally coupled scalar field system described 
by F{<j>)R coupling with F(<f>) = (l - £cp N ) R(N > 2) and the field potential, V{<£) = Vo(f>™. We 
use a generic set of dynamical variables to bring out new asymptotic regimes of the underlying 
dynamics. However, our dynamical variables miss the most important fixed point— the de Sitter 
solution. We make use of the original form of system of equations to investigate the issues related 
to this important solution. In particular, we show that the de-Sitter solution which is a dynamical 
attractor of the system lies in the region of negative effective gravitational constant Gjv thereby 
leading to a ghost dominated universe in future and a transient quintessence(phantom) phase with 
Gjv > around the present epocfQ. We also carry out comparison of the model with other competing 
^ models of dark energy such as galileon modified gravity and others. 

> 

I. INTRODUCTION 

OV 

\^ ■ Theories with a scalar field non-minimally coupled to gravity dubbed scalar tensor theories have been studied for 
• ' decades. The first well-known example of non-minimal coupling a la Brans-Dicke theory, was proposed in 1961 with 
an aim to match Mach principle with General Relativity [if. In this theory the gravitational constant is replaced by 
1 a scalar field <j> entering into the action in a specific combination with Ricmanian curvature as <f> 2 R. 

Followed by the Brans-Dicke proposal, other forms of scalar-tensor action were investigated, a well known example 
of a noii-minimally coupled system is provided by F(tfi)R coupling with F = 1 — £</> 2 . The cosmological dynamics of 
such a theory is rather rich and deserves attention. For a recent development in this direction, it is worth nothing 
■ that non-minimally coupled Higgs field due to a large coupling £ might give rise to a successful inflation Q which is 
* otherwise impossible. 

The non-minimally coupled scalar field system due to novel features are of great interest to dark energy model 
building pj-fTTj]. For instance, non- minimal coupling might allow phantom crossing and may give rise to cosmological 
scaling solutions of interest to models of dark energy. Phantom scaling solutions are generic features of a non-minimally 
coupled system with F = 1 — £ <fi 2 [l3| . 

In recent years, methods of dynamical system theory have been extensively used in cosmology for obtaining a 
general picture of dynamics for many cosmological models including those with a scalar field and modified gravity. 
The advantage of this method is in having some kind of " machinery " for deriving asymptotic solution using a simple 
programmed algorithm. This requires introduction of new set of variables in which the initial system can be rewritten 
as a system of first-order equations. On the other hand, there is a danger of losing some important solution as well as 
inability for this scheme to find transient regimes which can also be important. Nevertheless, a classification of stable 
asymptotic regimes, given by this method may be useful for understanding the underlying dynamics. The success and 



1 However, as demonstrated by Starobinsky in 1981, the ghost dominated universe, if exists, can not be accessed from the Universe we 
live in, we shall say more about this important result in the last section. 



2 



the limitation of the framework of dynamical systems applied to f(R) theory can be found in Rcfs.fTj] and fl5| |. 

We shall restrict our discussion to a polynomial functions F of the form F(<p) = 1 — £<fi N (N > 2) and power-law 
potentials for the scalar field, V (<$>) — Vo<t>" giving rise to generalization of models considered earlier (l3l fl5l - fT7| (see 



also Ref . [18| on the related theme) . The set of variables that we use can help in bringing out some generic features 
of the underlying dynamics and new asymptotic regimes missed in earlier studies. We should, however, note that the 
set of variables used in the present paper is not useful for the study of approximate Einstein regime in the system 
under consideration. For detailed description of this regime, other methods are required fl9j |. Secondly, our variables 
miss the existence of de Sitter solution and in order to investigate its existence and stability, we need to go back to 
initial variables to perform the analysis. 

In this paper, we investigate cosmological dynamics of non-minimally coupled scalar field system with specific 
functional forms of coupling and field potential using a convenient set of dynamical variables. We shall focus on the 
asymptotic regimes of the solutions of interest and reveal the important features associated with the de Sitter solution. 
We also carry out comparison of the model with other competing models of dark energy, in particular, the galileon 
system which is generically non-minimally coupled. 

II. EQUATIONS OF MOTION 

Let us consider the scalar field system with non-minimal coupling in the form, 
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m 2 Pl R-{g^^^ v + iRB^) + 2V{4>)) + S M , (1) 



where rn 2 pl = l/8irG = 1/k, and £ is the dimensionless parameter and Sm designates matter action. 
In a homogenous isotropic Fricdmann-Robertson- Walker universe with spatially flat metric, 

ds 2 = -dt 2 + a 2 (t)dl 2 , (2) 

the equations of motion which are obtained by varying the action ([I]) have the form, 

H " = I + + WHi'B'W + h2b ^)) + p) . (3) 

R = K (-<p 2 + 41/(0) + 3£(3£T05'(0) + ^B{4>) + 4> 2 B"{4>) + 4>B' (</>)) + p(l - 3w) J , (4) 



4> + 3H<P+^RB'(<j>) + V'(q)) =0. (5) 

where p and p are the energy density and pressure of the ordinary matter with p = ujp. 
From the standard form of equations 

1 ,f 
we read off the expression for Gn as the effective Newtonian gravitational constant 

° N = 87r(l-<£(0))' (6) 

and we shall use this definition of Gn hereafter. Ricci Scalar in chosen metric is given by, R = 6(2H 2 + H). For 
convenience, we shall use the system of unites with k = 6 . 

Dividing ©, 0) ,© by H 2 (l - 6£B(<f>)) and multiplying © by ££'(</>), we obtain, 

H 2 {1 - 6£B(») H 2 (l ~ 6£B((f>)) H(l - 6^5(0)) H 2 (l - 6£B(<f>)) ' 1 ' 
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R 6</> 2 241/(0) 54£<^B'(<£) 18£0 2 5"(c/>) 



B 2 B 2 (l -6£B(<j>)) H 2 {l-6£B((f>)) H(l - 6£B(<£)) H 2 (l - 6£B (</>)) 

lS&B'W 6 P (1 - 3aQ 

^2(1-6^(0)) H 2 (l- 6£B(0))' lJ 

WB'W 3&BW R eB' 2 (t) V'mB'W ( . 

H 2 (l - 6£B(<I>)) B(l -6£B (</>)) B 2 2(1 - 6£B(0)) H 2 (l- 6£B(0)) ' l J 

In the discussion to follow, we shall use the dimcnsionless variables, 

2 2V((f>) 6£4>B'(4>) 



ff 2 (l-6£B(<£))' V ff 2 (l-6£B(</>))' " ff(l - 6£B (<?!>) ) ' 

n = B 2 (1-6£B (({,))' (10) 
and the dimensionless parameters that depend on the specific form of functions £?(</>), V(<f>), 

(i-6^))' ' v{® ■ 1 ) 

Here ' denotes derivative with respect to </>. 

Using the new variables and parameters, we can now cast the equations of motion (J7J, (JSJ, (|9]) in the autonomous 
form, 

2i dx § _ _ _ _ tC _ , ^i' 7 " _ . 

x = 12X 2x( 2) +xz, 



B d\na z 6 

?/ tit/ , yz c Y 

dZ = * = *X+~-*&-*) + - 



i? dlna 6£A 6 

y4. d^L . , z , , , 

where wc have written equations for variables and parameters x, y, z, A, VL and dropped b & c as we would be 
interested in the cases of b,c = const (see Sec. IIIip . While writing equation for fi, we used the continuity equation, 
p + 3Bp(l + lu) = 0. Secondly, the higher derivative terms involving B and (f> in the equations of motion, 

x= Y= ± fl31 

- B 2 (l- 6^5(0))' ff 2 ' 1 J 

and Q get expressed through autonomous variables x,y, z and are given by, 



il = 1 — x — y — z, 

X( X ,y,z) = -f -T^^C-^ + ^+S + lz +3(1 -^-^-^ )(1 - 3 ^)' 

Y(x, y, z) = (^ 6x + U y + J_ (2b - ^) + 3(1 - x - y - z)(l - 3w)) . (14) 

In what follows, we shall investigate the autonomous system (jT2J for fixed points. We would specially be interested in 
stable solutions of interest to late time cosmic acceleration. For simplicity, we shall assume specific functional forms 
for functions B((f>) and V (</>). 
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III. STATIONARY POINTS AND THEIR STABILITY : B(4>) = 4> N , V(<j>) = V (j> n 

In this special case, b = B" (4>)<fi/ B' ((j>) = N — 1 and c = V'{4>)4>/V{(t)) = n. Our autonomous system was written 
keeping this simple case in mind. We shall find stationary points equating to zero the left-hand sides of the system (|12l) . 
Their stability will be established using the sign of the corresponding eigenvalues which we shall obtain numerically. 
We begin our discussion from the case, N ^ 2 (b ^ 1). In the case of N = 2, a simple additional relation exists 
between x and z, this case will be taken up separately in the Appendix. In the general case the form of this relation 
is more involved algebraically and its substitution into our system leads to cumbersome equations. We, therefore, 
prefer not to make use of it, we would rather check resulting solutions for consistency. 



A. The case of b 1(N ^ 2). 

After solving the system of algebraic equations obtained after equating the left hand sides of (fT2")l to zero, we find 
the following stationary points: 

1. Vacuum stationary line: x = 1, y = 0, z = 0, A € (— oo, +oo), fl = 

This is a vacuum solution for which the corresponding eigenvalues are given by, 

Ai = 3 - 3w, A 2 = 6, A 3 = 0, A 4 = 0. (15) 

As one of these eigenvalues is positive, the stationary line is unstable for any value of £ and w. 
The time dependence of scale factor can be found using the second relation in (fT3")) , 

Y stat = ^ = 6^2 + ^ 

Since Y sta t is constant, the above equation can easily be integrated to obtain the expression for a(t) 

i 

a(t) = a \t-t \ 2 a B ili , 

where t is integration constant. Since Y stat = — 6 for the stationary point under consideration, we finally have 

a(t) = a \t-t \K (16) 

We note that at all stationary points either t — > to or t — > oo. And all functions of time (a(t), <fi(t), H(t), p(t)) 
cither go to zero or infinity or become constants in this limit. For instance, the stationary line coordinate A has the 
following limiting values 



'0,if(f>(t) -> 0, 



N<f> ' 



tv , if (f> — (f)o — const. 



i-6^ J 

1). In case of 4>{t) — > and A — > 0, we find 4>(t) using the coordinate x of this stationary line. As (1 — 6^^) — > 1 for 
4>{t) — > 0, we have 

±4> = V^tIt H = ±V^7f , (17) 
0(t)=±-^gjL]n|£H 

6 

where 

n 6 

x s tat > 0, <j)' = const, t' — t' = y = const. 

a a 2 v* 5 tot 
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For the stationary line under consideration, x sta t = 1, Y sta t = — 6 which gives 

<Xt)=±±ln||^|, 



(18) 



Thus 4>{t) — > ±oo for t — > to or t — > oo. The obtained behavior of <p(t) contradicts the initial assumption that 
4>(t) — > and this case should therefore be discarded from the discussion. 



2). In case <j>(i) —> oo and A 
coordinates 



where 



-J|, we can find cf)(t) using the following combinations of stationary point 



■f,-f„i 

A stat 



i 2 2 

(pa 



H 2 N<p N N<P N a 2 
•,\2 



Istati 



fir = Nj st at{-) , 
where either -y sta t > 0, <j> N > or ~f sta t < 0, 



< 



(19) 



t'-t' = 



2 _ jalal 



const, j s tat = const, <f>' = const, . 



Since, for this stationary line, j sta t = x stat /A stat =1/A= -6£/iV, Y stat 

(2-JV)V6jf| 



-6, we have 



W)\' 



±- 



■ln 



t'-tn 



This solution exists for < i¥ < 2 
3). Finally, if 



oo = const, A = , again by using the combination 7 s t a t — { 

substituting it in the expression of 0(f) in eq. (|19[) . we get the same behavior, namely, 



± 



In 



t-*o 



(20) 



(21) 



and 



(22) 



Thus 4>{t) — > ±oo for t — > to or t — > oo and contradicts the initial assumption that <fi = <fio = const. 

= ~ 6f 



We therefore conclude that the stationary line reduces to a stationary point, A — 



We can now write the expression for the quantity Gm using the coordinates of the stationary point, 
A = B'(0)0/(1-6£B(0)) = N<t> N /(l-6$(f> N ) = N /((/>- N -6£), we find that 4> N = A/(N + 6(A) leading to 
G N = 6/8tt(1 - 6£ct> N ) = 6(N + 6£A)/8nN which is positive when £ > 0, A < -iV/6£ or £ < 0, A > -N/6£. 



Since the stationary line reduces to a stationary point for which A = — iV/6£ {4>{t) 
Gn — > 0. We also note that the vacuum solution does not contain parameters, w,n and N. 



oo), it follows that 



2. x = 0, y = 0,2 = 1, A 

The eigenvalues are, 



6+1 
6? 



Ai 



1 



N 



A, 



56 - c 



— , A 3 = 2-3u;, A 4 = 1. 



(23) 



6+1 N ' 6+1 
This point is unstable because A4 is positive for any £, to. We can compute a{t) by noting that Y sta t = 0, 

a(t) = a |t-io|* (24) 
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As for the expression for </>(t) , we use the following combination of autonomous variables at the stationary point 

Zatat _ jt_ _ a 



(25) 



where /? = = -gg and we have made use of the fact that (3 s t a t = const. For the vacuum stationary point under 
consideration, (3 s t a t = gff t = — = — ^ which finally gives 

0(t) = Mi-ior**. ( 26 ) 

Consistency analysis shows that this solution exists for iV > 2 and n < 5iV, so it never coexists with the solution 
discussed above under point 1. They both correspond to the situation in which the scalar field potential is negligible. 
Power indexes in this solutions do not depend on £, w, n, however the function <p(t) contains dependence on N. 
We also note that since A — — iV/6£ (<j)(t) — s- oo for t — > to) in the present case, Gn vanishes asymptotically at the 
stationary point. 

3. x = 0,y = 0,z = -l + 3w, A = = fl = 2 - 3cj 

The corresponding eigenvalues are, 

Ai 



(l-6)(l-3w) 



b + 1 

C (l-3o;) + 3(6+l)(l + a;) 

Aa = 6TI < ' 

for c> 3(6+ l),w < w < 1, 
2 

A 3 = -2 + 3w < 0, for u < -, 

A 4 = -1 + 3w < 0, for w < -. (27) 



where u; = - fj^g > | for 6 > 0, c> 3(6 + 1). 

As regions where A2 and A4 are negative do not intersect, this point is unstable (either a saddle or a repulsive node). 

1 — 3oj 1 — 3c 

6+1 ~~ N 



For obtaining a(t), <j){t), we note that Y sta t = 0, /3 = = i=*N. This tells us that, 



a(t) = a \t-t \\ (28) 

0(t) = ^oli-iol^. (29) 
The time dependence of the matter density we obtain using its continuity equation 

3(l + u) 

P (t) = po\t-to\ =po\t- to]- 31 '" (30) 
Power indexes of this solution do not contain parameters £, n and they depend on lj, N. 

It follows from the definition of Q that it must be positive to ensure that p > and Gn > 0. Since Q = 2 — 3w in 

the present case, we should have w < 2/3 to avoid a pathological situation. We also note that since A — —N/6£ 
(4>(t) — > 00), the effective Newtonian constant Gn vanishes at the stationary point. 



(6+l)(c(l-3w) + 3(6 + l)(w + l)) 

X = °' V = 2? ' 

_ 3(&+l)(^+l) b+1 N 

2c 2 - (b + l)(3(u; + l)(b + c + 1) + 4c) 
2c 2 
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The corresponding eigenvalues are, 
3(w + !)(&-!) 



Ai 

c 

3(c^ + 1)(6 + l) 2 + ic(u - 1)(6 + 1)) y/g + l)(/i(6, c)lu + / 2 (6, c) + 9^ 2 / 3 (6, c)) 
2 ' 3 4c(6 + l) 4c(6+l) 

A 4 = + D(1 + 6) > (32) 



where 



/i(6,c) = -210c 2 (6 + 1) + 162(6 3 + 1) + 19 26 2 c + 192c + 4866(6 + 1) + 3846c + 48c 3 , 
/ 2 (6,c) = 81(1 + 6 3 ) + (2436+ 17c 2 )(l + 6) + 174c(l + 6 2 ) - 16c 3 + 3486c, 
/ 3 (6,c) = (l + 6)(6+l + c)(9(6+l)-7c). 

and Ai, A4 are positive for — 1 < u> ^ 1 therefore this stationary point is unstable for any £. 

Analogous to previous points we get, 

_ 3(3(6 + l)(w + 1) + c(l - 3cj)) _ 3(3iV(w + 1) + n(l - 3w)) 

/? = 



then, 



c 




3(w + l) 3 


(w + 1) 


c 


n 


«(*) 


= ao|* 


0(i) 


= 0o|* 


P(t) 


= Pol* 



|*=. (33) 

We note that power indices of this solution depend on u>, N, n and do not depend on the coupling constant £. 

For 6 + 1 = c (N = n), the power index of functions a(t) and <p(t) is infinite and power-law solutions cease to exist. 
In this case the coordinates of the fixed point 4 are, x = 0, y = 2, z = 3(cj + 1), A = —^r: ^ = — (3w + 4). We note 
that < for u 6 [—1, 1] and that from the definition of i7 in (fT0|) . it follows that either p > 0, Gat < or p < 0, 
Gjv > which doesn't correspond to the real Universe. We also find in this case that Y stat = 12, /? = — 3 i ■ This 
tells us that H /H 2 = and <f>/(f> = /3Hq which is equivalent to 

a(t) = a e ffo( *- to) , 

ott tj. j. ^ 3 ^n( 1 +^)!*-*n) 

m = ^ H ^-^ = ^e — * , 

p(t) = poe -SH (i+u)(t-t )_ (34) 

As for the constants Hq and po, we substitute (|34|) in definitions (fT0|) of y, taking into account that N = n, y = 2 
and ft = —(3oj + 4). We then have 

2Vo<T 2y o 0o Ar e- 3H °( 1 +")( t -*°) V 

y 



n = 



H 2 (l-6&") Hq 2 A _ 6 ^ JV e _ 3 H (l +tl ,)(t- t0 )^ 3^ 

2/? 2/ 9oe- 3Ho(1+ " K *- to) po 



ff 2 (l-6£<^) ff() 2 A _ 6 ^ o iV e _3H (i +w )(t- t „)^ 3JT V£ 

(35) 



for t — > to when Hq > (or for t — > 00 when Hq < 0). 
Therefore, 



H 2 



Vq_ _Vq 
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where £ < 0, V > or £ > 0, Vb < 0, and 

po = -3r2iJo z 0o"£ 



Vo^o^w + 4) 



We next consider the behavior of quantities T = + T m and R — 6T e // = 8itGnT for the case N = n. Substituting 
the so obtained exponential solution, a(t) = aoe Ho ^~ to \ (f>(t) = <j>oc °' + «" ~ in (HJ, we find 

R 2 
Teff = — = 2if = consi, 

Q _ ^ (X e 3-ffo(l+«)(*-t ) 

8tt(1 - 6£</»o Ar e- 3ff «( 1 +")(*-*«)) 
T = + T m 

= e -3Ho(i-M(i-t„) ^ ^( 4Vo _ 27£(1 + w)# 2 + 27£(1 + w) 2 ff 2 ) + Po(l - 3w)) oc e -3ffo(i+")(t-to) (36) 
for t -> t Q , H > (or t -> oo, H < 0). 

It is therefore clear from the aforesaid that GW(i) grows as an exponent whereas T(t) decreases with the 
same rate thereby leading a constant product GnT. A remark about the exponentially expanding solution is in 
order. The solution though has features similar to de Sitter solution but does not really qualify for a true de Sitter 
as Gn is not constant in this case. We shall say more about this point in the discussion of the vacuum solution to follow. 

5. Vacuum solution 

This solution corresponds to the following fixed point, 

5 + 56-c 2(2 + 26-c) A 6+1 N _ „ 

x = 0, y=- , z= — ^ -, A = — = — -, n = 0. (37) 

y b+l + c b + l + c 6£ 6£' v ' 

The corresponding eigenvalues in this case are given by, 

2((6-l)(2(6 + l)-c)) ^ ^othj.^ 
Al = (6+l)(6+l + c) <0 > for &>1 and c< 2(6+1), 

A 2 = ^±^<0,for C <2(6+l), 
b + 1 + c 

5 _|_ f^fo Q 

A 3 = 1 < 0, for c < 5(6 + 1). 

(6 + 1)(3(6 + + 1) + c(7 + 3a,)) -2c 2 
4 (6+l)(6+l + c) 

for c < 2(6 + 1) when cj G [-1; 1], (38) 

The negativity of eigenvalues for this vacuum point show that it is stable for c < 2(6 + 1). Indeed, it is an attractive 
node for 6 > 1 and c < 2(6 + 1) for any numerical values of £ and u, see FigOTa). 
We note that, 

6c(5 + 56 - c) _ 6n(5iV - n) 
(6 + 1 + c)(6+ 1) ~~ (N + n)N ' 

2(2 + 26-c) _ 2(2N-n) 
(6+l + c)(6+l) ~ (N + n)N 

which gives rise to the following expressions for a(t) and 4>(t), 

a(t) = ao\t-t \^-^ K -^, 

<t>(t) = <j> \t-t \^ (39) 

This solution contains parameters TV, n but is independent of £ and tu. We note that power indexes in (f3"9")l arc 
negative for N < n < 2N and, therefore, a(t), <p(t) diverge leading to "Big Rip" singularity at t = to- This result 
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FIG. 1: The left panel (a) shows the Phase portrait for negative coupling constant £ = —1/2 for the fixed point 5, the 
fixed point is an attractive node. The right panel (b) displays the Phase portrait for positive coupling constant £ = 0.3961 
(N = 2, n = 7) for de sitter which is an attractive focus. 



2(;(2+„)-i) 



It — fn 1 2-11 obtained in the 



is generalization for the analogous vacuum solution a(t) = ao|t 
case of N = 2 [UBS. 

Let us further investigate the nature of the fixed point. For 6 + 1 = c (iV = n) and also for 2(6 + 1) = c (2JV = n) 
power indexes of functions a(t) and <fi(t) diverge and power-law solutions should transform into exponential ones. 
Indeed, the coordinates of the fixed point in these cases are given by, x = 0, y — 2, z — —1,A = — ||,f2 = and 

- J_ whereas 



1, z = 0, A = — g-j, f2 = respectively for N = n and 2N = n. In the first case, Y = 12, /3 = 



Y = 12, P = in the second. In case of N = n, we find that, jp- = 0, H = Hq = const and 



a e 



H (t-t ) 



„PH Q (t-t ) 



a(t) 

4>{t) = </) e^-' u ^ "■" = (p Q e 
We can find out H , using the definition of the coordinate y from (flU)) , 



(40) 



2^ nl 



2Fo0o JV e H °( t -*°) 



fl2 ( i_6e^) ^(i-e^o^ 



o(*-*o) 



/or f ->• oo , (i/ > 0) or t -> f (i? < 0). 
3iJo ^ 



(41) 



Thus, i? 2 = = where f < 0, F > or £ > 0, F < 0. 



In case of n = 2./V, the fixed point under consideration docs not lead to any physically admissible regime ( our 
numerical work shows that all trajectories in this case lead to oscillations near the minimally coupled field solution). 
It is interesting to note that in the absence of the standard curvature term in the action {F{<f)R = <fi N R), a family 
of de Sitter solutions exists (2(| for arbitrary values of field 4>- As for n = N, the solution corresponds to H = 
and w e ff = — 1 taken usually as the definition for de-Sitter solution tacitly assuming that Newtonian gravitational 
constant is a true constant of nature. In case, <f> is constant, the constancy of Gjv is trivially satisfied. However, in 
the model under consideration, we have an interesting vacuum solution with H = and an exponentially expanding 
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FIG. 2: This figure shows the real parts of eigen values for N = 2 and n = 5 versus £. The figure shows that the real part of 
two of the three eigen values is positive for £ < 0.1333 making the de Sitter solution unstable. Beyond the quoted value of £, 
real parts of all the eigen values remain negative. The two eigen values are complex conjugate of each other (their real parts 
coincide) in case of , 0.1333 < £ < 0.2068 and we see only two curves. Beyond £ = 0.2068, the eigen value curve splits into two. 



solution for <\> which corresponds to an exponentially decreasing/ increasing (depending upon the sign of Hq) the 
effective Newtonian gravitational constant, Gjv- The true de-Sitter corresponds to H = and (f> = const implying 
Gjv = const. Let us note that in this case, choosing O = 0, one might think to obtain dc Sitter solution but the 
latter implies H = 0. Clearly, this solution docs not qualify for a genuine dc Sitter. 

In fact, the true de-Sitter solution is not captured by the autonomous variables, we have used. In this case, the 
combination Ax + z 2 appearing in the denominator of autonomous system vanishes thereby telling us that autonomous 
variables used here are not suitable for the investigation of de-Sitter solution. We shall discuss this special case in a 
separate section to follow. 



IV. INVESTIGATION OF de-SITTER SOLUTION: ITS EXISTENCE AND STABILITY 



As mentioned earlier, our autonomous system is not suitable for the study of de-Sitter solution corresponding to, 
x = 0, y = 1, z = 0, A = — SI = as the combination, 4x + z 2 appearing in the denominator of dynamical system 
vanishes identically. In this special, we go back to the original variables and make use of the system of equations (|3|), 
Q, ([3]). Substituting, 7j = = = p = O and R = 12Hq 2 , for the de-Sitter solution in these equations we find, 



H Q 2 (1 - 6^ N ) = 2V fo n , Wq^NcPo"- 1 + IWo"" 1 = 0, (42) 

which gives important information for solution under consideration, 

2 V n(j) n - N N n 

H °- 6^' 00 -6Z(n-2NY (43) 

Clearly, this solution does not exist for n > 2N in general for an even n. We also note that for de-Sitter solution , the 

effective Newtonian gravitational constant, Gn = 6/87r(l — 6££J((/>)) = 3(27V — n)/8irN is positive only if n < 2N. 
The system of equations |(3J), Q, ([5) for p = which is of interest to us reduces to, 

H(l -6£B + 9£ 2 B' 2 ) = 3<p 2 (£B" -I) -3£B'(4:H<fi + V' + 6£B'H 2 ), 

4>{1-6£B + 9£ 2 B' 2 ) = -(3H<j> + V')(l -6£B) -^B'(-'i<p 2 + 12V + 9f + <p 2 B")). (44) 



In order to investigate the stability of the de-Sitter solution, we consider small perturbations \i & v around this 
background: H = Hq + [i and 4> = cf>Q + 1> m the dynamical system (|44|) . Consequently, we find the evolution equations 
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FIG. 3: The left panel (a) shows the evolution of (1 — 6£B((j>)) versus redshift z, where Gjv = 6/87r(l — 6£B(<j>)). The right 
panel (b) shows the evolution of w e j j versus redshift z. Both the figures correspond to the case of non-minimal coupling with 
N = 4, n = 9, £ = 1/2. We had chosen appropriate initial conditions, H = 2.034, <j> = -0.300,0 = -0.001 and £ = 0.5000 
(corresponding to w m — initially) to obtain w e ff ~ — 1 at the present epoch. The black dot in both the figures designates 
the present epoch which occurs in the regime of Gjv > 0, the effective Newtonian constant changes sign thereafter in future. 
The square marks the epoch where Gn turns negative. 



for perturbations 

M36£ 2 N 2 H Q fo 2N - 2 ) + /i(l - 6^00* + 9£ 2 iV 2 o 2Ar - 2 ) + u(3£VoNn(N +n- 2)q> N+n ~ 3 
+36£ 2 iV 2 (7V - l)F Vo 2iV - 3 ) + HmNHofo"- 1 ) = 0, (45) 

v(V Q n(n - l)0o"" 2 (l - 6^0^) - 6^ M^o W+n_a + 12£V N(N + n- l)^" -2 ) 

+P(3H (l - 6S<h N ) + 27eN 2 H Q( f> Q 2N - 2 ) + - 6£<h N + ^N 2 ^' 2 ) = 0. (46) 

In order to cast equations (|45p . (j4T>|) in a simple form, we introduce the following notations 

D = 36£ 2 iV 2 iJ o o 2JV - 2 , 

E = 1-6^0^ + 9£ 2 N 2 (j) 2N - 2 , 

F = 3ZV Nn(N + n-2)(j) Q N+n - 3 + 36£ 2 N 2 (N -l)H a 2 (j) Q 2N - 3 , 

G = l2^NH 4> N -\ 

Di - V n(n-l)fo n - 2 (l-6^o N )-^V Nncj) N+n - 2 + 12tVoN(N + n-l)fo N+n ~\ 

Eh. = 3i/ o (l-6^o JV ) + 27e 2 A r2 i/o0o 2Ar " 2 , 

Fx = 1-6^ N + 9£ 2 N 2 cj) 2N - 2 



(47) 



and make use of new variables s% = fi, S2 — s 3 = v. The system of equations (|45[) and (|46l) then acquires a simple 
form, 

Dsi + Efi + Fs 2 + Gs 3 = 0, (48) 
Dis 2 + E 1 s 3 + Fxi) = 0. (49) 
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In (1+z) 

FIG. 4: This figure shows the evolution of the equation of state parameter versus redshift z for a scalar field it in galileon 
model with Friedmann equation (|53 j) . The initial conditions were chosen carefully allowing us to be in the region of phase space 
where de Sitter is a stable fixed point [U]. The black dot on the curve designates the present epoch. 

Taking derivative of si, sa, S3 with respect to time we get the system of equations, 

d f g \ / \ 
~E ~e ~e\( s i \ 

1 \\ s 2 \ (50) 
-f -f J W 

which we have studied numerically. The numerical investigation of this system for b = N — 1 = 1,2,3,4,5,6 
and c = n = 1,2,3,4,5,6,7,8,9 shows that the de Sitter solution is stable when (we assume Vq = 1) 
2N + 1 < n = 5, 7, 9, —,(N — 2,4,6,...) and for positive values of the coupling, £ > £ > 0, where £0 de- 
pends on 6, c. We hereby quote few cases which show how the nature of the fixed points crucially depends upon 
numerical values of the coupling £: 

• For b=l (N=2) and c=n=5, the fixed point is 
saddle if < £ < £ , Co ^ 0.1333, 

attractive focus if £0 < £ < 0.2068, 
attractive node if 0.2069 < £ < 1; 

• For b=l (N=2) and c=n=7, the fixed point is 
saddle if < £ < £ , £0 ^ 0.0952, 

attractive focus if £0 < £ < 0.3999, 
attractive node if 0.4000 < £ < 1; 

• For b=l (N=2) and c=n=9, the fixed point is 
saddle if < £ < £ , £ - 0.0740, 

attractive focus if £0 < £ < 1; 

• For b=3(N=4) and c=n=9, the fixed point is 
saddle if < £ < £ , £0 ^ 0.0009, 

attractive focus if £0 < £ < 0.0014, 
attractive node if 0.0015 < £ < 1. 

The nature of cigen values of the perturbation matrix corresponding to eq. (|50p crucially depends upon the 
numerical values of b, c and £. For instance, the eigen values for the fourth point, (N,n,£) = 4,9,0.0012 are, 
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FIG. 5: In this figure, we plot the statefinder parameters for non-minimal model. The left panel (a) shows the time evolution 
of the statefinder pair {r, s} . The model converge to the fixed point (r = 1, s = 0).The right panel (b) shows the time evolution 
of the statefinder pair {r, q}, the model diverge at the point (r = l,q = 0.5) which corresponds to matter dominated universe 
and converge to the point (r = 1, q = —1) which corresponds to de Sitter expansion. The dark dots on the curves show values 
{ro, so} (left) and {ro, qo} at the moment when W e ff — — 1 from FIG. 3 (b) (right). Black squares show (r = 1, s — — 1, 0) (left) 
and (r = l,q= -1,0.5) (right). 



-47990.1409, -2286.0463 ± 1703.0817i which shows that de Sitter solution is an attractive focus. We numerically 
investigated the region of stability of the solution, see Fig|2] We did not find stable de Sitter solution in the range of 
n < 2N which has important implication for late time cosmology. 

The stable de Sitter solution exists in the region where the effective Newtonian constant Gjv is negative which 
means that graviton is ghost thereby leading to instability. Before approaching the attractor, the system passes 
through a phantom phase and parameters in the theory can easily be adjusted such that we obtain the observed 
value of equation of state parameter, w e ff ~ — 1 at present with Gjq > followed by a brief phantom phase before 
approaching the stable de Sitter fixed point ultimately pushing the ghost dominated regime to future (see Figj3J|. It 
is also possible to set the phantom phase at the present epoch. We have carefully managed to shift the ghost regime 
to future by adjusting the parameters in the model. We should, however, admit that such a model of transient dark 
energy suffers from ugly fine tuning problem. 



V. COMPARISON WITH THE COMPETING MODELS 



There are several popular models of dark energy which can comply well with popular observational data sets. 
The list of models is rather large, we shall restrict to commonly used models: ACDM, Quintessence with tracking 
behavior, Chaplygin gas, DGP and recently discussed galileon model which generically belongs to the class of non- 
minimal coupled models. 

One of the methods to distinguish between various models at the background level is provided by statefinder analysis. 
The Hubble parameter H and the deceleration parameter q constructed from the first and the second derivatives 
of the scale factor are often used to characterize a cosmological model. The statefinder method relies on the use of 
parameters constructed from third derivatives (and higher) of the scale factor [2^, [23[. For instance, in the lowest 
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FIG. 6: (color online) The left panel (a) shows the time evolution of the statefinder pair {r, s} for galileon (red) and Non-minimal 
(magenta) models. Both galileon and non-minimal models lie to the left of the ACDM fixed point (r = 1, s = 0). For both the 
models, s increases to zero from -1, whereas r first increases from unity to a maximum value, then decreases to unity. Both 
models converge to the fixed point (r = 1, s = 0) which corresponds to ACDM. The right panel (b) shows the time evolution 
of the statefinder pair {r, q} for galileon (red) and non-minimal (magenta) models. Both models diverge at the same point 
(r = l,q = 0.5) which corresponds to a matter dominated universe (SCDM) and converge to the same point (r = 1, q = — 1) 
which corresponds to the de Sitter expansion(dS).The dark dots on the curves show current values {ro,so} (left) and {ro,qo} 
(right) for different dark energy models. In all models, we have taken, Qom = 0.3 at the current epoch. 



order of this hierarchy, we have 



r - 1 



3(9-1/2) 



= 1 



1 w 
3wH : 



(51) 



which shows that (r, s) = (1, 0) & (1, 1) for ACDM and matter dominated universe respectively. In what follows, we 
shall briefly describe the aforementioned dark energy models followed by their statefinder analysis. 



VI. DARK ENERGY MODELS 



• Galileon model[2lj. 



Galileon is a massless scalar field, its Lagrangian apart from the kinetic term contains non-linear derivative 
terms responsible for switching off the modification to gravity locally. The theory is free from Ostrogradki 
ghosts despite the higher order derivative terms in the Lagrangian. The galileon action is of the form, 



(52) 



where {c,} are constants, /? is coupling constant of field with matter and L" [21| are Galileon Lagrangians. Li 
is linear in field and is often omitted assuming c\ — 0, L2 represents the standard kinetic terms, L3 = (c^tt) 2 !!]-^ 
is the famous Vainshtein term in which three it's participate; L4 and L5 are higher order Lagrangians. At least 
L4 is needed to obtain dark energy solution in the model(2l|. In this case, the evolution equations in a spatially 
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FIG. 7: (color online). The left panel (a) shows the time evolution of the statefinder pair {r, s} for Galileon (red), Nonminimal 
(magenta), DGP (blue), quintessence (Cyan) and Chaplygin gas (green) models. All models converge to the fixed point 
(r = 1, s = 0) which corresponds to ACDM.The right panel (b) shows the time evolution of the statefinder pair {r,q} for 
Galileon (red), nonminimal (magenta), DGP (blue) and chaplygin gas (green) models. All models diverge at the same point 
(r = 1, q — 0.5) which corresponds to a matter dominated universe (SCDM), and converge to the same point (r = 1, q = — 1) 
which corresponds to the de Sitter expansion(dS).The dark dots on the curves show current values {ro,so} (left) and {ro,(?o} 
(right) for different dark energy models. In all models, fiom = 0.3 at the current epoch. 



flat background have the form[2l| 



ZH 2 = p m + ^tt 2 - 3c 3 Htt 3 + yc 4 ff 2 7r 4 , 



2H + 3H 2 = -yiV 2 - c 3 n 2 if + -C 4 7T 3 



(3H 2 it 



2Htt + 8Hn , 



-e 2 (3Hn + tt) + 3c 3 tt (3ff 2 rr + Htt + 2H%) - 18 Ci HTT 2 (3H 2 7r + 2Hi\ + 3Hn 



The energy density, pressure for the scalar field tt and the Friedmann equation are, 



Pn 



H 2 



2 -tP - 3c 3 Hn* + yC 4 iP7r 4 , 



-7T- + C 3 7T^7f - -C 4 7T J + 2HtT + 8HlT 

8irG r , .o -, Pit 
—Z-[P0m{l + Z) +p 7r \] Wir = 



As demonstrated in Rcf.[21(, there is a region in phase space which necessarily contains a stable de 
solution. In FigJH we have plotted equation of state parameter (with favorable initial conditions) 
shows that the latter ultimately approaches its de Sitter value. 

• Dark energy with a constant equation of stateJBJ. 

The Hubble parameter is given by, 

H{z) = H Q [fi 0m (l + zf + n DE (l + z) 3 ^) 1 1 



(53) 
(54) 

(55) 
(56) 

(57) 

Sitter 
which 



(58) 
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where (ACDM) corresponds to w = — 1 and for Quiessence [22| w = constant ^ — 1. 
• Non-minimally coupled scalar fleld[13|. The energy density and pressure in this case are given by, 

= ^ 2 + V( ( f>)+3t;(Hj ) B'(4>) + H 2 B(0)), 

P = \ft-V{<t>)-Z (2H^B / ^) + <j> 2 B' / (4 ) ) + ^B'(4 ) ) + (2H + 'm 2 )B(4 ) ) 



The Hubble parameter for this model has the form, 

8irG 



H 1 



[/?0m(l + zf + P<t> ] 



P4> 

W<f, = — 



(59) 
(60) 

(61) 



Quintessence^, Q- 

In this case, there is a variety of models which can lead to desired cosmic evolution with a carefully chosen 
field potential. All the quintessence models can be classified into two groups: Thawing and Freezing or tracker 
models. In case of thawing, cosmic evolution crucially depends upon initial conditions where as the tracker 
dynamics is free from the choice of initial conditions and we shall focus on the latter. Bearing that in mind, we 
have the expression for Hubble parameter, 



tt2 87rG r 



^) 3 + 



V{cj>) oc (f>~ a ,a > 1 



(62) 



Chaplygin gas [291 ). 



H(z) = H 



Oom(l + zf + 



^Orn A 



+ {l + zf 



1/2 



k V B 

• The Dvali, Gabadadze , Porrati (DGP) model |3 

H(z) 



K= POm/vB, A — B < K 



2 / 1 — ^Qn 



H 



l-f2 0m \ (1 — fioir 

+ WOom(l + Z) A + 



(63) 



(64) 



In FigJSJ we have shown the time evolution of the statefinder pairs {r, s} and {r, q} for non-minimally coupled system. 
The model converge to the fixed point (r = l,s = 0) which corresponds to ACDM. The right panel (b) shows the 
time evolution of the statefinder pair {r, q}, the model diverge at the point (r = l,q = 0.5) which corresponds to a 
matter dominated universe (SCDM), and converge to the point (r = l,q = —1) which corresponds to the de Sitter 
expansion (dS). 

For comparision with other models, we have displayed our model along with other competing models in Figs.((6])& 
(0. Since, to the best of our knowledge, statefinder analysis for galileon is not done elsewhere, we have in Fig.© 
displayed (r, s) trajectories for non-minimally coupled scalar field system and galileon field which is generically non- 
minimally coupled system. FigJT] displays statefinder trajectories for non-minimally coupled scalar field along with 
other popular models of dark energy including the galileon field. The figure shows that all models diverge from the 
same point (r = l,q = 0.5) which corresponds to a matter dominated universe (SCDM) and converge to the point 
(r = l.q = —1) which represents de Sitter univcrse(dS). Figj6] shows how clearly the two non-minimally coupled 
systems get distinguished on the (r,s) & (r,q) planes. In case of galileon, we focussed on the sector of the model 
which contains stable de Sitter fixed points. FigJT] allows to distinguish these two models from other popular models 
of dark energy. It may be noted that we have re-scaled appropriately the statefinder variables of the model under 
consideration in order to strike comparison with other models. It is interesting to note that the behavior of statefinder 
pairs for galileon and Non-minimally coupled scalar field systems resembles with that of the braneworld models studied 
in Rcf . |32| . 



VII. CONCLUSION 



In this paper, we have revisited cosmological dynamics of non-minimally coupled scalar field system. We presented 
detailed investigation of dynamics of the underlying system in case of F((f>)R = (l — £B(<fr))R coupling with B(<f>) = <fi N 
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and V{<j>) = Vocff 1 (N > 2) using a convenient set of autonomous variables. We studied asymptotic regimes of solutions 
in the model. In case of the vacuum solution, we found a very interesting solution for which H = {w e tt = —1) 
and the scalar field <f> increasing exponentially giving rise to exponentially decreasing Gn — effective Newtonian 
gravitational constant. Such a solution does not qualify for de Sitter for which Gn should be held constant. 
The autonomous variables used in this paper though convenient in general but miss certain important features of the 
dynamics. The description fails to capture the de Sitter solution as the combination of variables Ax + z 2 appearing in 
the denominator of the autonomous system vanishes identically in this case. The investigation of this solution took us 
back to the original variable in the evolution equations. We found that in case of de Sitter solution, Gn > provided 
that n < 2N. On the other hand, our numerical investigations showed that the solution under consideration is stable 
only for n > 2N + 1 (we checked for lower values of N > 2). For initial conditions of matter dominated universe, 
the system enters the phase of acceleration consistent with observation at present followed by a brief phantom phase 
thereafter which continues till de Sitter is reached, see Fig|3£b) (which is an attractive focus of the dynamics, see 
FigJIJb)). During the phantom phase Gn changes sign from positive to negative thereby making the universe ghost 
dominated in future, see Figj3] It is possible to set parameters in the model such that phantom phase occurs at 
present epoch corresponding to Gn > compilable with observed values of the equation of state and fractional 
density parameter pushing the ghost dominated phase to future which no body has yet seen. Incidentally, similar 
features of equation of state are shared by the braneworld model discussed in Ref . (33j . 

An important remark about the ghost dominated universe with Gn < is in order. It was demonstrated by Starobinsky 
in 1981 that any attempt of crossing over to the regime with Gn < makes the Friedmann universe unstable at the 
transition point. Any tiny perturbation over homogeneity and isotropy grows large at the point of transition. Taking 
into account the quantum effects in the scalar sector does not effect this result 1 . The ghost dominated universe with 
Gn < 0, if exists, will be separated from our real universe by a region of large curvature of the order of Planck value. 
Thus it is impossible to cross the point where the effective gravitational constant changes sign due to the formation of 
generic anisotropic curvature singularity at this voint\3A^ . 

The de Sitter solution for the case with N = 2 also shares the aforementioned features. In other cases, we have 
shown that the solutions obtained earlier are continued for values of N > 2. We have investigated in detail the 
asymptotic regimes of these solutions in all cases including the one corresponding to N = 2 

We have shown that the non-minimally coupled scalar field system can account for late time cosmic acceleration. 
We should, however, emphasize that dark energy in this scenario appears as a transient phenomenon which involves 
extra fine tuning. 

Last but not least, we used statefinder analysis to compare the model with other rival models of dark energy, in 
particular, with galileon model which generically a non-minimally coupled system. The model under consideration 
is clearly distinguished on the (r, s) and (r, q) planes from other competing models of dark energy. We focussed our 
attention on sectors which contain stable de Sitter solution in case of the model under consideration as well as the 
galileon modified theory of gravity. 

A final remark about the ghost dominated evolution in the model is in order. In the framework of the simple set 
up of non- minimal scheme discussed here, there exists no consistent de Sitter solution such that Gn remains positive 
throughout the evolution. It would really be interesting to explore generic functional forms of the coupling function 
and the field potential to check for a well behaved de Sitter solution. 
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Though the higher order curvature corrections to Einstein-Hilbert action motivated by Starobinsky in different context could help to 
avoid the catastrophe and land the theory ghost free. 
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VIII. APPENDIX 

A. The case of b = 1(N = 2). 

As mentioned before, this case was considered earlier in other references. However, discussion of asymptotic regimes 
was missing and important issues related to the structure of de Sitter solution were not discussed. Thus we have 
retained this subsection in the appendix which is also useful for comparison of our results with b = l(N = 2) case. 

As mentioned above, there exists a simple relation between x and z for N = 2. Indeed we have, 



H 2 {l-^B{4>)) ff2(l-6£B(0))(6fB'(0)) a #l-6£B(0)) 72? A 
Now we substitute (p5|) and b = 1 in (fH|) , 

z 2 

= 1 7TT — V — Z, 

72?A y 
z_ _ 1 fz 2 
2 ~ 1 + 18£ 2 A V12' 

+3e 2 ^(i - ^4 " 

1 / z 2 

1 + 18? A \12?A' 



. 2 



(2 - 3c£) + 3(1 - - y - z){l - 3w) ) . 



Excluding from (TP2"]) x, f2, A, Y - using (j6"6")l . we obtain following system for b = 1 (A = 2) 

'(WTWA) {wa^ - 1} + 6y{2 - 3c?) + 3(1 - Wa - y - ^ - ^) - 2 ) 



6£,4 



(65) 



-3, - - 1) + ytMA + c) + 3^(1 - ^ - „ - ,)(! - 3.))) 

■£a-* Ui + Is? a) (Wa & ' V + 6 ^ 2 - 3 ^ + ^ - Wa ~ y ~ ~ 3 ^) ~ 2 ) + 



c 2 . 



A' = ^ + Az. (67) 

We enumerate the stationary points and corresponding solutions of (|6"T)) , 
1. de-Sitter solution 

y = l,z = 0,A = -^,n = 0, 

Y = 12,0 = 0, 

a(i) = a e ffo(t - to) , 

<£(t) = 0o, (68) 



where 0o = ±./ fiff "_ 4 -, , H = ±\/ — V ° n f^ — are obtained from the system ([3])- (O for H = (f> = O = p = 0. 



6J(n-4) ' u — V 12 5 

This is de Sitter solution which exists also for A > 2. We have not solved the autonomous system in this case as 
the combination 4x + z 2 appearing in the denominator vanishes in case of de Sitter solution with coordinates x = 0, 
y=l,z = 0,A = -^,Q = 0. 
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y = 0,z = 12^ + 2^36^ -6^, 

Y = 6(12^-1 + 2^/6^-1)), 
/3 = -6£ - V^-l), 

1 

a(t) = dolt-tol 3-12 *- 2 ^ 6 * -1 ', 

6C + V6t(6{-1)" 

</>(*) = ^ojt — *o| 3 - 12 «-V6«(^-D ; (69) 



y = 0,z = 12£- 2^36^ -6£, 



F = 6(12£-1-2V6£(6£-1)), 
P = -6f + A /6£(6£-l), 



-6;+V6{(s;-i) 

= ^olt-^ol 3 - 12 ^ 2 ^^, (70) 

Solutions corresponding to stationary points 2, 3 were found in fl6| . where the cosmological model with non- minimal 
coupling F(<f>) = ^4> 2 and the effective potential V(<f>) = \(j) n (without Einstein term in the Lagrangian) was 
investigated. 



4. 



?y - * 4e+c-i !j4 " 

n = 1 2g(l - 3a.) 2 4g(l - 3tj) 
3(4£ + o;-l) 2 4£ + o,-l' 
= 3(l-3o,)(o,-l) _ 2£(1 - 3oQ 
4^ + w- l ' P 4£ + u> - 1 ' 

2(4{ + ^-l) 

a(t) = a |t-t | 3 ^+ 16 «- 3 , 

4g(l-3m) 

0(t) = 0oI*-ioK+ lf *- 3 , (71) 



This solution also was found in [3] . 
5. 



(6£o;(c - 2) - 2£(c + 6) + 3 - 3w 2 ) 
y ~ 2c^ ' 

_ 6(1 + 0,) 1 

3(1 - f (c + 2))w + c 2 £ + 3 - 7c£ - 6£ 



3(6 + c-3o,(c-2)) _ 3(1 + 0,) 
c c 
a(t) = a |i-io| 3( ™- 22 " (1+ " ) , 
4>{t) = <f> \t-t Q \^, 

p(t) = polt-tol*^- (72) 
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We note that the analogous solution exists in the case N > 2 (see (|33p ). It also was obtained in [l6j . For n = 2 
the power index of functions a(t) and </>(i) becomes infinity and power-law takes the exponential form. In this case 
the coordinates of the point 5 arc y = 16 ^ g 3 + 3 " 2 ; z = 3( w + = -^,fi = 3(1 ~ 4 ° 4 ^ 16C+3 and also, Y = 12, 
j8 = — 3 ( 1 +"^ ) . This means that, = 0, H = Ho = const, % = [3Hq giving rise to the following 



„H (t-t ) 



a(t) = aoe 
0(t) - 0oe^ o(t - 4o) = o e " , 

p(t) = i0oe -3i?o(i+")(t-*o). (73) 
In order to find the constants Hq and po, we substitute (|73|) in definition (TTUJ) of coordinates y, ft taking into account 



that n = 2, y = lb ^ 3 + 3 " and n - 



16£-3+3^ 2 anA n _ 3(l-4g)^-164+3 

2Y O 0™ 21/ (/>o 2 e- 3ifo(1+ " )(t - to) Y 



y = 



ff 2 (l-6£0 2 ) H o 2 (l-6C0n 2 e- 3ff o(i+-)(*-to)) 3 # 2 £ 

2p _ 2 A) e- 3g °( 1 +"X t - t °) po 

' i? 2 (l-6£0 2 ) ff 2 (1 - 6^ 2 e- 3ff o(i+-)(t-to)) 3i/o 2 0o 2 C 



for t ^ to when i?o > 0(or for i — > oo when Ho < 0). 
Therefore, 



#n 2 = 



3£y 3(16£- 3 + 3w 2 ) ' 



where £ < 3( ^ } , Y > or f > 3(1 r"' } , Vq < and 



16 ' u ■ ^ S 16 

po = — 3f2i2o 0o £ 



2 ^ 2 , 2W> 2 (3(l-4£V-16£ + 3) 



6. 



16$ - 3 + 3w 2 

(6c 2 £ 2 - c 2 £ + 8c£ - 48£ 2 c - 120£ 2 + 56g - 6) 
6(£(2 + c)-l) 2 

2 ^-,a = -±-« = o, 



£(2 + c)-l' 3$' 
y _ 3(10c$ - c 2 g - 4) £(c-4) 



£(2 + c)-l £(2 + c)-Y 

2({(2+n)-l) 

a(t) = a |t-t |«"- 2 )("-^, 

0(t) = 0o|i-to|^ (75) 

This solution in the limit £ — > oo gives (|39p [l6j . 

For n = 2 and n = 4 power indexes of functions a(i) and <f>(t) become infinity and power-law solutions transform into 

exponential ones. We find that coordinates of point 6 in these cases are, y = 96 ^ 4 ~ 3 ^ j" 3 , z = — 4 ^£ 1 , A = — Q = 
and y = 1, z = 0,A = -i, ft = accordingly for n = 2 and n = 4. Then Y = 12, /3 = - for n = 2 and Y = 12, 
(3 = - for n = 4. For n = 2 we obtain -§2 = 0, H = Ho = const and 

o(t) = a e Ho{t - t0 \ 

4>{t) = 0oe 0H ^-^ = foe^^ . (76) 
Using the definition of the coordinate y from (Tl7J|) analogous to previous point we find the constant Hq (n = 2, 

y = 9 m-1$ 3 ) and hcncc H ° 2 = ~m = - w$-~m +3 y where e e (-°°. °) U(i &)» v > o or $ G (o, ±) u(&)> 

Vo < 0. For n = 4 wc get | = and 

a{t) = a e H " {t " to) = ao, 

0(f) = O = 0, (77) 
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where values Hq = 0, (f>o = were found by substitution a(t) = ciQe 110 ^ *% (f>(t) — (f>o, n — 4 in the system of 
equations J3])-([5]). 
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